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Abstract 

The time evolution of a system where a uniform and classical SU(3) color electric field and 
quantum fields of quarks are dynamically coupled with each other is studied focusing on non- 
perturbative pair creation and its back reaction. We characterize the color direction of the electric 
field in a gauge invariant way, and investigate its dependence. Momentum distributions of created 
quarks show plasma oscillation as well as quantum effects such as the Pauli blocking and interfer- 
ence. Pressure of the system is also calculated, and we show that pair creation moderates degree of 
anisotropy of pressure. Furthermore, enhancement of pair creation and induction of chiral charge 
under a color magnetic field which is parallel to the electric field are discussed. 

1 Introduction 

Study of non-perturbative pair creation from a classical electric field, which is known as the Schwinger 
mechanism [Ij, has a long history and wide range of applications (see Ref. [2j for a recent review). One 
of those applications can be found in studies of relativist ic heavy-ion collisions, where the Schwinger 
mechanism has been used as a mechanism of matter formation from a color flux tube [3J. The color 
flux-tube model assumes that a strong color electric field is formed in a beam direction just after two 
nuclei collide and pass through each other ^ . Formation of longitudinal color electric fields is also 
predicted in the framework of color glass condensate Therefore, particle production due to the 

Schwinger mechanism attracts renewed interest ^ [9l [TOl [HI [T2J . 

Under these circumstances, getting an understanding of how an initial electric field and created 
particles evolve in time is of prime importance. To properly describe the time evolution, calculating 
vacuum persistence probability or pair creation probability, which were first derived by Schwinger, 
is not sufficient jl3], and an electric field should be treated as a dynamical variable rather than a 
background field controlled by hand, i.e. back reaction should be taken into account. There have been 
considerable numbers of studies treating back reaction; the ones based on a kinetic theory [3l [HI ITS] 
and the others on quantum field theory [131 [161 [TTl [HI [19] . To our knowledge, however, field theoretical 
treatment of the back reaction problem under a color electric field has been lacking. Therefore, in 
this paper we investigate the pair creation of quarks under a color electric field incorporating back 
reaction. 

In studies of physics under non-Abelian electromagnetic fields, SU(2) theory has been often used 
for simplicity. In the case of SU(3), however, a new feature arises: anisotropy in color space. It has 
been shown that an SU(3) color electric field has two independent directions and it is characterized 
by two gauge invariant parameters: one of them is determined by its field strength and the other 
is related with the color direction of the field [20^, More generally, an SU(A/'c) color vector has 
[Nc — l)-independent directions in color space, and physical contents can generally depend on a color 
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direction of an electric field [22j. In this paper, we deal with SU(3) color electric fields and examine 
the color direction dependence. 

Not only new features which arise in non-Abelian fields, we also analyze phenomena whose essence 
is common to the Abelian case. Collective motion of created particles which couples to an electric 
field shows plasma oscillation. During this evolution, several phenomena are observed: suppression of 
pair creation or annihilation of the particles due to the Pauli blocking, damping of the electric field, 
and rapid oscillations in the momentum distribution of the created particles due to interference. We 
shall give an analysis of these phenomena to advance an understanding of physics in pair creation. 

We take a uniform color electric field as an initial state. Pressure of this initial state is quite 
anisotropic: the longitudinal pressure is negative and the transverse pressure is positive. Therefore, if 
local thermalization is achieved starting from the flux-tube initial condition, isotropization of pressure 
should be needed during the time evolution. However, the full understanding of a thermalization 
process in heavy-ion collisions has not been obtained. In this paper, we examine the role of pair 
creation for the isotropization of pressure as a first step to understand a mechanism of thermalization 
in heavy-ion collisions. 

One of remarkable differences of the color flux tube given by the color glass condensate from that 
in the original flux-tube model is the existence of a longitudinal color magnetic field in addition to 
an electric field [?]• It has been shown that a longitudinal magnetic field enhances pair creation of 
fermions and speeds up the decay of an electric field in the previous paper [13j. We extend it to the 
quark pair creation under a longitudinal color electric and magnetic field. 

Furthermore, we study induction of chiral charge due to pair creation under a magnetic field. Since 
the chiral anomaly is a semi-classical effect where the quantum aspect of a gauge field is unnecessary, 
we can also apply our framework to study the chiral anomaly due to pair creation. The relation 
between pair creation and the chiral anomaly has been also studied in Refs. [241 125] . Emergence of 
a nonzero chirality in heavy-ion collisions attracts interest in the context of the chiral magnetic effect 

m- 

The remainder of this paper is organized as follows. In the next section, we shall explain the 
Abelianization of a color electromagnetic field, and introduce the parameter characterizing the color 
direction of the field. Although this formalism is essentially the same as that given in Ref. [26], we 
make the existence of color direction dependence clearer with the help of the method in Refs. [20^ 2^. 
In Section|3j we introduce time-dependent particle picture to describe the time evolution of the system. 
Then, we shall show our numerical results in Section [Ij Time evolution of momentum distribution 
functions of created quarks, color current density, electric field strength and pressure of the system are 
displayed and discussed. Color direction dependence of the results is also examined there. In Section 
[5} effects of a longitudinal magnetic field, i.e. enhancement of pair creation and induction of chiral 
charge, are discussed. 

2 General framework 

Quark pair creation incorporated with back reaction is described by the following Lagrangian density 

C = i;{irD^-m)i;-^F^,F-'''', (1) 

where ^ is a quark field and color indices i (z = 1, 2, • • • , Nc) are omitted. We assume for simplicity 
that each flavor has the same mass m, and flavor indices are also omitted. The number of flavor is set 
to be Nf = 3 throughout this paper. 

The covariant derivative and the field strengths are defined in terms of a background gauge field 
Af, as 

D, = d^ + igT''Al, (2) 
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= d,Al - d.Al - gr'^A^Al, (3) 

where is the generator in the fundamental representation of gauge group SU(A/'c), and f^^^ is the 
anti-symmetric structure constant (a, 6, c = 1, 2, • • • , N'^ — 1). 
The equations of motion for and '0 now read 

D^F-^^ = ^^7"T>, (4) 
{i-f^D^ - m) ^ = 0. (5) 

Because we treat the gauge field as a classical background field, the charge current operator 
^'07^r^'0 in Eq.Q is replaced by its expectation value g{i(;^^T^tl;) in the following. These coupled 
equations govern the pair creation and its back reaction. 

We restrict the background field to spatially homogeneous and Abelian-like one, which is expressed 

as 

= F^,n«, (6) 

where Fj^j^ is an Abelianized field strength and is independent of the space coordinates, is a constant 
vector indicating a color orientation of the electromagnetic field, and is normalized so that n^n^ = 1. 
This field strength is given by the gauge field 

The relation between F^^^ and is just the same as that in the Abelian case: 

F^u = d^A, - d,A^. (8) 

Under the Abelianized gauge field Q, the covariant derivative ^ reads — igrf'T^A^. 
Because {n^T^)- is an Nc x Nc hermit ian matrix, we can diagonalize it by a unitary transforma- 
tioifj global gauge transformation): Un°'T°'U'^ — diag(tyi, • • • ,wn^). Then, the equations of motion 
(|4j) and ^ can be rewritten as 

Nc 

^,F^^ = J2'^^9W^^^P'). (9) 

1=1 

[i-f^{d^ + iwigA^)-m]r^O (i = l,2,--- ,iV,). (10) 

These equations have the same form with those in the Abelian theory, except the presence of the 
weight vector Wi. Hence, we can solve the coupled equations Q and (10) as in the Abelian case [ 13] . 
Each quark field couples with the background field via a coupling constant Wig. 
In the case of SU(3), a diagonalized n^T^ may be expressed as 

Un^'T^'U^ = T^cos9-T^sm9, (11) 

because = diag(l/2, -1/2, 0) and = diag(l/2\/3, l/2\/3, -l/\/3) are diagonal. (We represent 
as = where is the Gell-Mann matrix.) Then, Wi are expressed in terms of 9 as follows 

1 71" 1 StT 1 371" 

" 71 COs(e + -), W2^^ COs{e + W3^^ COs{e + y ). (12) 



*^Diagonalization itself is always possible by local gauge transformation even if a color electric field is not spatially 
uniform. In that case, however, Wi acquire space dependence, and the method used in this paper is no longer applicable. 
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Figure 1: A graphical representation of Eqs.(12). The triangle of dotted lines corresponds to = 



and is an usual weight diagram for a fundamental representation of SU(3). Each Wi is given by the x 
coordinate of each vertex of the ^-rotated triangle. 



The angle 6 is related with the Casimir (gauge) invariant [201 [2T]. our parametrization*^ , the 
relation is 

sin^S^^SCs, (13) 

where C2 = [d^^^n^n^n^]^ is the second Casimir invariant for SU(3) and characterizes a direction of 
the color electromagnetic field in a gauge invariant way. 

Fig{l]is a graphical representation of Eqs.(12). Each Wi is given by the x coordinate of each vertex 
of the rotated triangle (weight diagram). Owing to the symmetry of the diagram, it is sufficient to 
take the angle 9 restricted to < < |. 

3 Canonical quantization in background fields 

To describe the time evolution of the system, we introduce an instantaneous particle picture and 
quantize the field using that particle picture. Of course, a definition of particle in the presence of 
a pair-creating background is rather ambiguous. Thus, the instantaneous particle picture should be 
regarded as our working definition to describe the system evolution. To calculate field quantities such 
as current and energy density, we need to identify and subtract the contribution from the Dirac sea 
at each time. We can do this just by the normal ordering in terms of the instantaneous particle basis, 
and then energy conservation is automatically guaranteed [T3]. Therefore, one can interpret defining 
a particle picture at each time as a means to properly regularize field quantities, and can consider 
particle number at intermediate time as a byproduct. Nevertheless, we will show time evolution of a 
particle number or momentum distribution function defined through the instantaneous particle picture 
in the following, because it behaves in a physically reasonable manner, and helps us understand the 
dynamics of this system. 

First, we consider only the electric fields. We will introduce the magnetic fields in Section [5j We 
set an initial state with no electric field and switch it on at time t = 0. An initial field strength is £"0 
and its direction is along the z-axis. The gauge = is chosen so that 

^(t) = o (t<0), 



dt 



= -Eo, (14) 
t=o 



^The relation between in this paper and those in Refs. |20lET| is ^ = 0Ref.[2T] + f = f^Ref. 20 + 
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and A\ — A2 — 0. After switching on, A'^{i) evolves according to Eqs.Q and (10). 
A quantized quark field may be expanded as 



3=^,1-^ 



i inf 



(15) 



where ttp^^ [bp^] is the annihilation operator of a particle [antiparticle] with momentum p and spin-5 

satisfying the anti-commutation relation {ay ^j^^f} — {^p,s5^qT/^} — ^ij^ss'^^{p ~ q)? and ±^l^ps{x) 
are classical solutions of the Dirac equation ([To|. The superscript 'in' distinguishes the initial condition 
for ±tl^ps(x): at t < 0, -i-i^ps{x) [-i^ps{x)] is identical to the positive [negative] energy solution of the 
free Dirac equation. We set the state to be in- vacuum |0,in), where no particle exists initially and 
which is defined by a^'^^\0, in) = ^'j^^^lO, in) = 0. 

At t > 0, ±^l^ps(x) evolve under influence of the electric field and become superposition of a positive 
and negative energy (frequency) state. To describe the pair creation dynamically, we introduce a 
time-dependent particle picture by decomposing the field operator into positive and negative 

frequency instantaneously: 



(16) 



where +'0p5^°^(x) [_'0p5^°^(x)] is a positive [negative] frequency solution of the Dirac equation under 
the pure gauge = As{t = to). Instantaneous particle picture is defined by ap^^(t) and b'^p^si'^)- Of 
course, ap^^(t) and 6p^s(t) agree with ttp ^ and b^^^ at t < 0, respectively. The particle picture at time 
t and that of the in-state are related by the time-dependent Bogoliubov transformation: 



i inf 



of which coefficients satisfy |Q;p5(t)P + |/3ps(^)P — 1 and are given by 



(17) 



<s{t)SHP - q + WigAsit)) = /d3x+Vpi*^^(x)+C;(x), 



x+#^(x)-V'^?(x). 



(18) 



A quark pair distribution function is defined by 



f;,{t) = (0,in|a;,t(f)a;,,(t)|0,in)^ = (0, in|6!_Tp,(t)6Lp,(t)|0, in) 



(27r) 



(27r) 



V 



V 



(19) 



where V is the volume of the system. With the help of Eqs.(17), we can rewrite fps{t) in terms of the 
Bogoliubov coefficients: 

fUi) = Ksm'- (20) 

The expectation of the charge current operator, which is regularized by the normal ordering is 



^'^W = ^O'^^l • ^i^S^' : |0,in 

i=l,2,3 

d^P Pz rz 



i=l,2,3.=U ^ ' ^ i=l,2,3«=t,i 



(21) 
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where the transverse mass mx = ^jm? and the anomalous distribution 



|0,in) 



(270 



Re [e-2-^*a*p,(t)/3^,(t)] (22) 



are introduced *^ The first term of the last expression of Eq.(21) is a conduction current, which is 
caused by movement of real particles and the second is a polarization current, which is generated by 
variation of electric dipole density. 



4 Pair creation under electric fields 
4.1 Plasma oscillation and Pauli blocking 



We have solved the coupled equation ^ and (10) numerically with the help of Eqs.(20), (21) and (22). 



The results are shown in Figs J2| [3] and [4] The time evolution of longitudinal momentum distributions 
with fixed transverse momentum is exhibited in Fig J2] and that of transverse momentum distributions 
with fixed longitudinal momentum is in FigJS} Fig|4] shows the time evolution of color current density, 
electric field strength and quark number density. The parameters are set to be a = = 0.01 or 
a = 0, and ^ = 1,0 = 0. To concentrate our attention at first on general features which are common to 
the Abelian case, only the distributions of "blue" quarks, whose effective coupling to the background 
is when = 0, are presented. Hereafter, all figures are shown in the dimension-less unit scaled by 

After the switch-on of the electric field, particles are created with approximately longitudinal 
momenta. Their occupation number is approximately exp (^—Trm^/lwilgEo^, which accords with the 
form expected from the semi-classical tunneling calculation [27], and the transverse momentum dis- 



tribution exhibits a Gaussian-like form * In particular, particles in zero mode (m = 0,pt = 0) are 
created as many as possible under the restriction of Pauli's exclusion principal fp^si'^) ^ 1? so that 
their occupation is exactly equal to one. Therefore, the longitudinal momentum distribution with 
a = and = [Figj2]^b)] shows squarish form. 



Oscillating factor e~'^^^p^ in the anomalous distribution is canceled out in total because both aps(t) and I3ps{t) 
contain the factor e*^^*. 

*^ This is not exactly Gaussian. An exact Gaussian distribution is realized in the case of a constant electric field (no 
back reaction). Under a constant electric field, we can calculate the distribution analytically and can show that the 
distribution approaches to Gaussian in an asymptotic region where pz ^ oo. 
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(a) « = 2^ = 0.01 





(b) a = 



Figure 3: Transverse momentum distribution of "blue" quarks, g = l^p^/ y/gEo = 1.0,^ = 



After created, particles [anti-particles] are accelerated to the direction of WigEi [—WigEi]. They 
generate charge current in the direction of n^, which decreases the electric field according to the 
"Maxwell's" equation Q. Then, the direction of the electric field is flipped at some time and the 
particles start to be decelerated. Repeats of this process result in oscillating behavior of the longi- 
tudinal momentum distribution, the charge current and the electric field, which is known as plasma 
oscillation. 

Other than plasma oscillation, which is a classical dynamics, also the Pauli blocking, which is 
a quantum effect, plays a role in the time evolution of the longitudinal momentum distributions. 
Because of Pauli's exclusion principal, a particle blocks the subsequent pair creation at the point 
where it locates in phase space. In particular, if a value of a distribution exceeds 1/2, not only pair 
creation is suppressed but also pair annihilation occurs (see the second term of the right hand side of 
Eq.([32|)). Therefore, the distributions get dented when particles cross the line of = in momentum 
space. (Notice that pair creation and annihilation can happen only in the vicinity of the line of 
Pz = because we now take only a classical gauge field.) The effect of the Pauli blocking is the most 
conspicuous in distributions with mx = 0. In this case, because the occupation of the created particles 
is equal to one, particles are totally annihilated when they cross the line of = [See Figj2]^b)]. 

Collecting the facts above, the momentum distribution function can be approximated by the 
following equation: 

fi,{t) exp [-J^^) [p,{-WigA,-p,)) . (23) 

This equation roughly replicates the distributions obtained by the numerical calculation. Using this 
empirical and analytic expression of the distribution, we analyze the numerical results and study the 



parameter dependence of current and particle number density. Equation (23) is exactly correct for the 
zero mode with — 0. In contrast, it loses its accuracy for higher mx modes. However, it does not 
matter for our purpose because pair creation of those modes is strongly suppressed compared with low 



mx modes. Although Eq.(23) is neither an exact one nor obtained by some systematic approximation, 
it reproduces the numerical results with sufficient accuracy for a rough analysis and is useful because 
of its simple structure. 



Using Eq.(23), the momentum integrations in the particle number density and the conduction 



current can be done and we obtain 



N{t) 



jTT^Eo \Mt)\ (24) 



jr\t)c^-^Eo l^idfe ^oAsit). (25) 

^ ^'^ 1=1,2,3 
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Figure 4: Time evolution of the color current density, the electric field and quark number density. 
5 = 1,^ = 0. 



In the derivation of Eq.(25), the approximation Pz/^p — sgn(^^) {m? <C qEq) has been used, which 



makes Eq.(25) overestimate. Equation (25) can explain the plasma oscillation obtained by the nu- 



merical calculation. Neglecting the polarization current and using Eq.(25), we can solve the Maxwell 
equation = jz and obtain the oscillating electric field 



where its frequency is 



A3(t) = - ^sinVtt 
Ez{t) = Eocosnt, 



(26) 
(27) 



(2vr) 



mgre 



(28) 



Let us introduce the time tc when the electric field first reduces its strength to zero. This tc gives a 
typical time scale of variation of the electric field. In the present approximation, tc is given as follows 



_ 1 27r _ TT 



(27r) 



(29) 



This gives a smaller value of tc than the result of the numerical calculations because the current (25) 



is an overestimate. However, their discrepancy is less than factor 2 and Eq.(29) correctly describes 



the parameter dependence of tc obtained by the numerical calculations. Equation (29) tells us that 



the stronger the electric field or the larger the coupling constant, the shorter time the electric field 
takes to vanish. 



Substituting Eq.(26) into Eq.(24) leads the oscillating behavior of N{t) while the actual N{t) 
[FigQc)] does not show oscillation. This discrepancy is because of the incorrectness of the approxi- 
mation ([23]) at t > tc. For explanation, suppose particles having a positive charge Wig > 0. They get 
positive momentum at first by acceleration of the electric field, and after tc they go into the negative 
momentum area in the momentum space because the direction of the field is flipped at t = tc. The 



modeled expression (23) does not describe particles plunging into the negative momentum area. If 
the particles are those in the zero mode mx = 0, they totally disappear when they cross the line of 



Pz = due to the Pauli blocking and thus Eq.(23) is correct at all time. If mx 7^ 0, however, they 



do not totally disappear and the simple approximation (23) becomes incorrect after tc- That is why 



Eq.(24) is not correct at t > and the actual N{t) shows saturation rather than oscillation. Under a 



strong magnetic field, however, N{t) does oscillate in time [FigfT3[c)] and the present approximation 
becomes correct because only the lowest mx mode mainly contributes (See Section [s]). 
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Figure 5: Comparison between the polarization and conduction current, a = 0.01,^ = 1,0 = 0. The 
electric field is also plotted for reference. 



4.2 Polarization current and damping of the electric field 



The current density (21) consists of two parts: the conduction current and the polarization current. 
As discussed in the previous subsection, the conduction current is associated with collective motion 
of particles, i.e. plasma oscillation. In contrast, the polarization current is related with microscopic 
processes of pair creation. It is induced at an instant when pair creation happens. In a classical view, 
when a particle pair is created and becomes on-shell, a distance between them is nonzero except the 
case that their transverse mass is zero. Therefore, the pair creation process can be interpreted as 
creation of an electric dipole, so that it generates the polarization current. (See Ref. [13] for more 
detailed argument.) 

Because of its origin, the polarization current is expected to be induced in the same direction with 
the electric field. In other words, Ez and jz^^ have the same sign. Therefore, the polarization current 
always reduces the electric field strength through the Maxwell equation = —jz- That is, the 
polarization current causes damping of the electric field. 

Nevertheless, the polarization current obtained by our numerical calculation shows irregular be- 
havior [Figjs] rather than simple one expected from above argument based on a classical view. This is 
because the system undergoes a complicated evolution due to plasma oscillation and the Pauli block- 
ing, and furthermore the anomalous distribution, which is included in the integrand of the polarization 



current (21), is sensitive to phases of the Bogoliubov coefficients. What is especially remarkable in 
the plot of the polarization current is a peak just after switching the electric field on. Because the 
electric field is turned on suddenly, the vacuum experiences prompt polarization. That is why the 
large polarization current is induced at an early time. 

Although the naive expectation based on a classical view does not hold in this dynamic system, 
damping of electric fields actually happen. However, it is slight and hard to be recognized in Figjij^b) 
because the polarization current is far smaller than the conduction current (except the time just after 
the field is switched on) in our parametrization. If behavior in longer time is calculated, damping 
of the field would be observed. Indeed, damping behavior can be recognized in the result under a 
magnetic field which is parallel to an electric field [Fig|T2|^b)], because a magnetic field speeds up the 
time evolution of the system (see Section [5]). Besides it, the damping would be quickened by strong 
coupling. 

Before closing this subsection, let us emphasize that a zero transverse mass mode does not con- 
tribute to polarization current. It is evident from the explicit expression of a polarization current, 
which contains the factor mx (the second term of Eq.([2T|)). In a classical view, the distance between 
a massless pair is zero when they are created, so that its electric dipole moment is also zero. That is 
why no polarization current is induced by zero modes. 
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Figure 6: Time slices of the longitudinal momentum distribution in Figj2]^a). 
4.3 Interference 

One may notice that the longitudinal distribution shows rapid oscillations after tc [Figj2]. Fig J6] shows 
time slices of the longitudinal distribution with a — 0.01,^ = 1 and = 0, whose whole picture is 
Figj2]^a). Before tc 15/ in this case), the distribution is smooth. After tc, particles created 
before tc start to cross the line of Pz — 0. Then, the distribution begins to show rapid oscillations. 
These oscillation have been observed also in earlier works where back reaction is taken into account 
[T6l[l7]. These can be interpreted as interference between particles created before tc and those created 
after tc- 

Because the problem with back reaction cannot be treated analytically, we deal with pair creation 
under a constant electric field to explain the interference noted above. The problem without back 
reaction can be solved analytically and an explicit expression for distributions is available [13j. Under 
a constant electric field, however, particles are accelerated to one direction, so that there is no event 
that two distributions overlap in momentum space nor interfere. Therefore, we suppose an initial state 
with a distribution /o(p) instead of a vacuum with no particle |0,in). (In this subsection, we omit 
spin and color indices for simplicity.) This initial distribution would join a distribution of particles 
created from an electric field, and they would interfere. That initial state is expressed by a two-mode 
squeezed state as 

i/o)=n-^p"'^'e^p 

p 

where AAp and F(p) are related with the initial distribution /o(p) as follows 

|i-(p)P = r^' ^p = r47;^ ^^^^ 

l-/o(p) l-/o(p) 

Notice that a phase of F(p) is irrelevant to /o(p)- Because now /o(p) is given by hand, there is no 
criterion to decide a phase of i^(p). However, if /o(p) is a distribution of particles created from an 
electric field, a phase of F(p) is automatically determined by its time history of evolution. 




V 



|0,m), 



(30) 
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(a) neglecting the interference term 




(c) F(p) 



/o(p) „ivlleEo 
1-/0 (P)^ 



(b) F(p) 




(d) F(p) 



/o(p) ^-wl/eE^ 
l-/o(p)^ 



Figure 7: Distribution functions under a constant electric field with the initial distribution /o(p) 
ie-(p/V^+5)2^ a = 0.01,e = l,pT = 0. 



Letting this system evolve under the electric field which is given by the gauge — (0, A(t)), we 
obtain the distribution function at time t\ 



U{t) = {h\aUt)a^{t)\h) 



-pvv-pv-yuu/ y 

= /o (P + WigA{t)) + {1 - 2/0 (p + WigA{t))} /p(t) (32) 
+ 2 {1 - /o (p + WigA{t))}Re [ap(t)/3;(t)F(p + WigA{t))] . 

The first term of the right hand side represents the initial particles accelerated by the field. fp{t) in 
the second term is the distribution of particles created from the field, which is defined by Eq.([l9]). The 
factor {1 — 2/o (p + WigA{t))} expresses the effect of the Pauli blocking: initial particles suppress the 
subsequent pair creation, and pair annihilation occurs if the initial occupation exceeds 1/2. The third 



term describes interference between the initial particles and those created from the field, *^ Now ap{t) 
and /3p(t) are the Bogoliubov coefficients giving the distribution fp{t). Although the distributions 
fp{t) and /o(p) are independent of the phases of the Bogoliubov coefficients or F(p), the third term 
is sensitive to those phases. Because of this term, the distribution shows rapid oscillations when the 
two distributions overlap in momentum space. 



^This interference term has been missed in the previous paper [13j 
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As an illustration, we plot in Fig j?] the distribution (32) with the initial distribution /o(p) = 
under a constant electric field which is switched on at t = 0. For simplicity, results 
in Abelian theory, which is obtained by replacing Wig with e in ([32]), are shown. The parameters are 
chosen as a = 0.01, e = 1. For reference, we show in Figj7|(a) the result in which the interference 
term is neglected. We can see clearly the effect of the Pauli blocking. In the presence of the third 
terra, the distributions show oscillations of which contour is along classical trajectories of particles 
[Fig|7[b),(d)]. What is remarkable is Fig^c), in which a phase factor of F{p) is set to be e^Vzl^^^ ^ in 
this case, oscillations seen in FigjTj^b) disappears. That is to say, the factor e^vil^^^ cancels the phases 
of Q:p and /S*. This is reasonable because the phase factors of and /3p are expected to be e^*^^^ and 
^-iSci ^ respectively, in which is the classical action of a particle under the constant electric field: 

fe^oo). (33) 

Hence, if one want to reproduce the situation in which a distribution of particles created from the 
electric field interfere with a distribution of particles which are also created from the electric field, 
it is natural to assume F(p) has the same phase factor with «p/3p, i.e. ^-^Vzl^^^ ^ A distribution in 
such case is plotted in FigjTj^d). Rapid oscillations similar to those seen in the distribution with back 
reaction are obtained. 

Note that like as the polarization current, the interference term has no contribution from a zero 
transverse mass mode. That is because the occupation np^s{t) = \l3ps{t)\'^ of a zero mode takes the 
maximum value 1, and there is the constraint lo^p^ P + |/3psP = 1, so that Q;p(t)/3*(t) is for this mode. 

These rapid oscillations are a key ingredient for seeming irreversibility of time evolution [28]. 
Because of this interference as well as a polarization current, effective dissipation of energy from the 
electric field to the quantum fields happens. For example, the number density in FigjlFc) shows 
saturation, and there is seemingly an arrow of time, although the equations of motion Q and ^ 
hold time-reversal symmetry. In contrast, if there is no polarization current nor the interference term, 
effective dissipation would not happen and evolution of the system would be periodic in time. We shall 



demonstrate it by studying pair creation under a strong magnetic field in Section 5.1 (see FiglTSfc)). 



4.4 Color direction dependence 

As shown in Section [2| direction of a color electric field in color space can be characterized by the 
Casimir invariant C2 or the angle in a gauge invariant way. Therefore, physical values may generally 



depend on C2 or 0. Casimir dependence of transverse distributior^^ under a constant color electric 
field has been studied in Ref. [29j. However, as far as the author knows, color direction dependence 
incorporated with back reaction has not been investigated. 

Figjs] shows the color direction dependence of the longitudinal momentum distributions (a = 
0.01,^ = 1,Pt — 0). We can understand the behavior of these distributions based on the values of 
the effective coupling Wig. For example, when = 0, the distributions of "blue" and "red" quarks are 
symmetric in momentum space because wi = —W2, and "green" quarks are not at all created since 
ws = 0. Concerning the momentum distributions, the result greatly depends on the color direction 
of the electric field. This is a matter of course because the distributions of unconfined quarks are not 
color singlet object. 

Color direction dependence of the charge current density, the field strength and the quark number 
density is shown in FigJOj Unlike the momentum distributions, the color direction dependence is 
very small for these quantities. This is because the current and the particle number are obtained by 



Their definition of transverse distribution is different from us. It actually must be called pair creation probability 
rather than transverse distribution. 
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Figure 8: Color direction dependence of the longitudinal momentum distributions, a — 0.01,^ = 
1,PT 0. 



summing up the colors, so that they depend on 9 only through the elementary symmetric polynomials 
of Wi. The relations between 9 and the elementary symmetric polynomials of wi are 



2 I 2 I 2 



W1W2WS 



sin 3^. 



(34) 
(35) 

(36) 



The first equation is the trace of Eq.( 11 ), the second is the trace of the square of Eq.( 11 ) and the third 
is the determinant of Eq.( 11 ). Because wi+W2-\- and wl-\-W2+ are independent of 9 and only 



W1W2WS has the dependence on 9, the field quantities depend on the color direction parameter 9 very 
weakly. 

As an illustration, see Eqs.(24) and ([25|). Because our attention is now on the regime where pair 



creation strongly happens, i.e. rm? / gE^ <C 1, the factor exp ^— i^^j^^^ ^ is nearly equals to 1 and its 
^-dependence is negligible. Therefore, leading dependence on 9 comes from the factor (wig)'^ or 



However, the factor {wigY brings no ^-dependence because of Eq.(35). That is why color direction 
dependence of the number density is very small. In the case of the current density, although the 
factor \'^ig\^ does bring ^-dependence, one can show that the ^-dependence is numerically small 
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Figure 9: Color direction dependence of the current density, the electric field and the quark number 
density, a = 0.01,^ = 1. 



by substituting Eq.(12) explicitly into \uJig\^- It varies between 0.24^^ and 0.25^'^ 



4.5 Pressure 

The initial state with the longitudinal electric field = (0, 0, E)n^ is quite anisotropic: longitudinal 
pressure is Pl = —E'^/2 and transverse pressure is Pt = £'^/2. In contrast, pressure is locally 
isotropic in a locally thermalized quark-gluon plasma phase. Hence, if longitudinal color electric fields 
are formed in the initial stage of heavy-ion collisions and if local thermalization is realized in the 
later stage, isotropization of pressure must be achieved during the evolution of the system. In this 
subsection, we examine how the initial anisotropic pressure evolves under the influences of pair creation 
and its back reaction. 

Pressure generated by quarks is given by the expectation of the symmetric energy-momentum 
tensor: 



(0,in| : : |0,in) = --(0,in| : ^ (7^ 9 ^ + 7^ d^jtp: |0,in) 



(37) 



= diag(f,PT,PT,PL), 

which is regularized by the normal ordering. Notice that (O^^) is diagonal because we are in the center 



of mass frame. Inserting Eq.(16) into Eq.(37) and using Eq.(17), one can express £ (energy density; 
Fig, 10), Pl and Pt in terms of the momentum distribution fps{t) and the anomalous distribution 



Pl = 2Nf 



1 — 1,2,3 s — '1^,4' 



i=l,2,3 s=^,i 

Pt -2NfY. E 

2=1,2,3 s=^,i 



UJr, 



1 



(27r)3 2 



(t) 



(38) 
(39) 
(40) 



The terms containing fps{t) are contributions from real particles, which are common to classical theory. 
What is peculiar to our quantum field theoretical calculation is the presence of the terms including 
Qpsit), which exhibit pressure generated by pair creation processes. 



14 



10 20 30 40 50 



Figure 10: Energy balance. 
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Figure 11: Time evolution of pressure, a 
transverse, respectively. 
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(b) total pressure 



0.01,^ = 1,^ = 0. L and T denote longitudinal and 



One can confirm that the relation 



(0,m|:e^„:|0,in) 



£-2Pt 
2Nf 

i=l,2,3 s=tA 

m(0, in| : : |0, in) 



(2^ 



(41) 



holds 

The result of our numerical calculation is shown in FigfTT 
and 9 



0. 



The parameters are set to be a = 0, ^ = 1 
Due to pair creation and subsequent acceleration by the electric field, the particles 



generate positive pressure in the longitudinal direction, whereas transverse pressure by the particles 
is relatively small because the particles are accelerated only to the longitudinal direction. As pressure 
by the particles increases, pressure of the electric field is weakened owing to back reaction. The 
negative longitudinal pressure at the initial time is compensated by pressure generated by the particles. 
Although full isotropization is not achieved because this system is collision-less plasma, degree of 
anisotropy is moderated due to pair creation. 



5 Effects of magnetic fields 

In the color glass condensate framework, it has been shown that longitudinal color magnetic fields are 
produced in addition to color electric fields just after heavy-ion collisions [SI [7]. In this section, we 
study the effects of a longitudinal color magnetic field on pair creation. Let us consider the situation 
such that both electric and magnetic fields are along the z-axis in configuration space and along 
in color space. This setup realizes a gauge field configuration with nonzero topological charge 
pa^pa^^ _^ Q ig predicted by the color glass condensate. 

*^Trace anomaly is not included in our treatment because we do not consider quantum gauge fields. 
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Because this system is spatially uniform, induced current is also uniform, so that the magnetic field 
is not changed by back reaction. Therefore, the magnetic field acts as only a catalyst in this situation. 
To treat back reaction to a magnetic field, a non-uniform field must be considered. However, it is 
beyond the scope of our present paper. 



5.1 Enhancement of pair creation 



Under a longitudinal magnetic field, pair creation is enhanced, so that the time evolution of the system 
becomes faster due to the enhanced back reaction [13j. This phenomenon is caused by appearance of 
Landau levels and spin-magnetic field interaction. 

Due to interaction with the magnetic field, transverse momentum of created particles is discretized 
into Landau levels, and degeneracy between modes with spin parallel and antiparallel to the magnetic 
field, which are denoted by t and ^ each in the following, is broken: 



where 



n 

n+ 1 



— > 2ns\wi\gB, 
(n = 0,l,---) 



for s =t 
for s 



(42) 



(43) 



The expressions of field quantities under the longitudinal magnetic field are available from those under 

+ 2ns\wi\gB 



no magnetic field by the replacements (42) and 



UJ- 



(PpT 



\wi\gB 
2tt 



(44) 



n=0 



where p in ujp^s is the abbreviated expression for {pz^ Notice that the transverse mass with = 
is independent of while those with higher depend on B. Therefore, a strong magnetic field 
makes particles in higher modes "heavy" and suppresses their pair creation. In contrast, creation of 
{us — 0)-particles is not at all suppressed. 

Not only the creation of the lowest mode is not suppressed, the magnetic field enhances field quan- 
tities such as current and total particle number. That is because the number of modes degenerating 
in a unit transverse area is proportional to B. 



For example, the charge current (21) is replaced by 



j,{t) = 2Nf J2 E ^^^^ 

i=l,2,3 s=t,^ 



27T 



27r 



(45) 



In a strong magnetic field gB > gE ^ m^, this may be approximated as follows 

Pz 



jz{t)c^2N^ 



f J2 

=1,2,3 



\m\9B f dp^ 
2ti / 2ti 



{t) + 



m 



-9. 



p^,n=0,t 



(46) 



which is proportional to B. 

Fig{T2] presents the magnetic field dependence of the current, the electric field and the particle 
number density (a = 0.01,^ — 1^9 — 0). The color current density and the quark number density are 
indeed enhanced by the magnetic field, and as a consequence the frequency of the plasma oscillation 
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Figure 12: Magnetic field dependence, a = 0.01,^ = 1, 6^ = 0. 
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Figure 13: Mass dependence under the magnetic field. B = E^^g = 1, = 0. 
increases. We can estimate tc, which is a typical time scale of back reaction, in the same way as that 



in Section 4.1, and obtain 



7T 

tc= - 

2 



(2jr) 



tanhvr-^ 



B 



(47) 



This is a decreasing function of B. 

Because only the lowest mx mode mainly contributes under a strong magnetic fielc *^ , the results 
in the massless case (a = 0) and those with nonzero mass (a = 0.01) show drastic difference [Fig|l3 
whereas in the absence of a magnetic field, there is no considerable difference between them [Figj4 
The particle number density is especially remarkable: Its time evolution becomes nearly periodic in the 
massless case. Similar behavior of a particle number has been obtained by Iwazaki [25] in a different 
way. Emergence of this periodic behavior is consistent with our contention that the polarization 
current and the interference of distributions are the source of apparent irreversibility of time, because 



as noted in Sections 4.2 and 4.3, the zero mode contributes to neither the polarization current nor the 
interference term. 



5.2 Chiral anomaly 

The magnetic field has another effect: induction of chiral charge. Because the lowest Landau level 
(n = 0) is occupied only by particles with spin parallel to the magnetic field (s ^t)? the balance 
between spin-^ and ^ is broken and thus chiral charge is induced. 

We consider the chiral charge density given by the following equation 



Qs = (0,in| : ^7075^ : |0,in). 



(48) 



*^ However, contribution from higher modes is not completely zero if magnetic field strength is finite. In particular, 
the peak seen in the current density just after switching on the field does not disappear even in the case of massless 
particle under the magnetic field, although this is caused by nonzero uit modes. 
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Figure 14: Time evolution of the chiral charge, g — 1^ B — E^^O — {). 



In the chiral limit (m = 0), the chiral charge equals the difference between the number of right-handed 
particles and that of left-handed particles. 

Substituting Eq.(16) into Eq.(48) and subtracting a divergent part by the normal ordering, we can 
express the chiral charge density in terms of the distribution functions: 



^ WjgB r 



i=l,2,3 



dpz 

27V 



Pz 



m 



:0,t(^) 



(49) 



Note that only the lowest mode contributes to Q5. Contribution from the higher modes is totally 
canceled out because they are occupied with both the spin-t and ^ modes. When m 7^ 0, the anomalous 
distribution also contributes to the chiral charge. 

The time evolution of the chiral charge obeys the Adler-Bell-Jackiw anomaly equation ^30l [31]. 
The expectation of the anomaly equation extended to the QCD case now reade 



*9 



j^Q^{t) = ^E{t)B + 2mQ^{t\ 



(50) 



where the term 2mQ^ describes explicit breaking of the chiral symmetry due to the mass term, and 
Q5(t) = (0,in| : ^A^s^A : |0,in) 



(51) 



=1,2,3 



In Fig, 14 , each term of the time integral of Eq. ( 50 ) is plotted. Equation ( 50 ) is indeed satisfied. Due 



to plasma oscillation, the chiral charge density shows oscillation. As noted in the previous subsection, 
time evolution is nearly periodic in the chiral limit [Figfl4|^a)]. This curve of Q5 becomes exactly a 
sine one in the limit of a strong magnetic field B ^ Eq [25]. In contrast, if m 7^ 0, the term Q5 grows, 
and the periodicity seen in Q5 is destroyed. 



6 Summary and discussion 

We have studied the time evolution of the system where quantum quark fields and a classical color 
electric field are dynamically coupled with each other. By the procedure of Abelianization, we have 
reduced the equations of motion to the same form as in the Abelian case. Effective coupling strength 
of each quark with the Abelianized electromagnetic field depends on the direction of the field in color 
space. Therefore, the longitudinal momentum distributions of created quarks greatly depend on the 
color direction of the gauge field [Figjsj. However, we have found that the color current density and 

*^ The divergence of the chiral current is dropped because of the space homogeneity. 
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total particle density are rather insensitive to the color orientation of the field since those quantities 
are obtained after summing up all the color components [Figjo]. 

We have obtained the oscillating behavior of the longitudinal momentum distributions, and esti- 



mated its typical time scale tc [Eq.(29)] using the empirical analytic expression for the quark distri- 
butions. This tc also gives a time scale of net particle production. The increase of the quark number 
density is concentrated at initial time and gets saturated after tc [Figjij^c)]. Let us roughly estimate tc 
with the parameters motivated by the color glass condensate; m 0, A^/ = 3 and qEq 2^ ^ IGeV^, 



where Qg is a saturation scale. Then Eq.(29) yields 

tc - 3/g fm. (52) 

Note in passing that the dependence of tc is also small, 2% effect, at most, in the present case. 
This tc is not drastically small compared with the typical time scale of the initial stage of heavy-ion 
collisions, such as thermalization time expected from hydrodynamic calculations < 1 fm (e.g. [32]). 
However, this estimation incorporates only the effect of quark pair creation. If gluon pair creation is 
taken into account, tc would become much smaller. Pair creation of gluons is stronger than that of 
quarks because (i) gluons are boson, so that they are Bose-enhanced rather than Pauli-blocked, and 



(ii) effective coupling of gluons to an electric field is larger than that of quarks However, applying 
the method in the present paper to gluon pair creation, we encounter a difficulty of infrared divergence. 
Furthermore, under a magnetic field, one particle energy of a gluon in the lowest Landau level becomes 
imaginary and instability occurs, which is known as the Nielsen-Olesen instability [33l I34 j . Thus, we 
need to resolve these difficulties in order to treat both the quark pair and gluon pair creations in our 
framework in a unified manner. 

In addition to the plasma oscillation, we have pointed out that the response of vacuum to the 
electric field involves several phenomena such as the Pauli blocking, damping of the electric field 
and rapid oscillations in the quark distribution. In particular, it has been revealed that the cause 
of the rapid oscillations is the interference between the particles. Because of this interference, the 
distribution becomes sensitive to the phase of the Bogoliubov coefficients. Furthermore, we have 
found that apparent irreversibility of time evolution is induced by nonzero transverse mass modes 
through the polarization current and the interference. 

Also the time evolution of pressure of the system has been calculated. We have shown that 
the initial anisotropy in pressure is moderated by pair creation and back reaction even though our 
treatment is a mean field one. If effects of collisions among particles are taken into account, in other 
words, if quantum gauge fields are introduced, the system would be more isotropized. It is hoped to 
investigate these effects, as well as effects of non-uniformity of electromagnetic fields. 

Finally, we have discussed the effects of a magnetic field on quark pair creation. We have shown 
that a magnetic field enhances pair creation because of the emergence of the Landau level and the 
spin-magnetic field interaction. As a result, the longitudinal magnetic field speeds up the decay of the 
electric field. This mechanism may have significance in the context of heavy-ion collisions. 
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